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Abstract 

In  this  paper  we  analyze  the  beam  combining  and  atmospheric  propagation  of  high-power  lasers 
for  directed-energy  (DE)  applications.  The  large  linewidths  inherent  in  high-power  fiber,  and  to  a 
lesser  extent,  slab  lasers  cause  random  phase  and  intensity  fluctuations  occurring  on  sub¬ 
nanosecond  time  scales.  To  coherently  combine  these  high-power  lasers  would  involve 
instruments  capable  of  precise  phase  control  and  operating  at  rates  greater  than  -'10  GHz.  To  the 
best  of  our  knowledge,  this  technology  does  not  currently  exist.  This  presents  a  challenging 
problem  when  attempting  to  phase-lock  high-power  lasers,  which  is  not  encountered  when  phase¬ 
locking  low-power  lasers,  for  example  mW  power  levels.  Regardless,  we  demonstrate  that  even  if 
instruments  are  developed  that  can  precisely  control  the  phase  of  high-power  lasers;  coherent 
combining  is  problematic  for  DE  applications.  The  dephasing  effects  of  atmospheric  turbulence 
typically  encountered  in  DE  applications  will  degrade  the  coherent  properties  of  the  beam  before 
it  reaches  the  target.  Through  simulations,  we  find  that  coherent  beam  combining  in  moderate 
turbulence  and  multi -km  propagation  distances  has  little  advantage  over  incoherent  combining. 
Additionally,  in  strong  turbulence  and  multi-km  propagation  ranges,  we  find  nearly 
indistinguishable  intensity  profiles  and  virtually  no  difference  in  the  energy  on  the  target  between 
coherently  and  incoherently  combined  laser  beams.  Consequently,  we  find  that  coherent  beam 
combining  at  the  transmitter  plane  is  ineffective  under  typical  atmospheric  conditions. 
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I.  Introduction 

Advances  in  solid  state  lasers,  especially  slab  and  fiber  lasers,  have  made  them  candidates  for 
directed-energy  applications  [1,2].  To  achieve  the  power  levels  needed  for  these  applications  it 
is  necessary  to  combine  a  large  number  of  lasers  into  a  single  beam.  The  combined  beam  is  then 
propagated  many  kilometers  through  a  turbulent  atmosphere.  Recently,  much  effort  has  been 
expended  in  coherent  beam  combining  [2-7].  Coherent  beam  combining  refers  to  matching  the 
phases  of  multiple  lasers  either  in  the  transmitter  plane  [2,3]  or  the  target  plane  [4].  However, 
there  are  a  number  of  important  issues  to  be  considered  before  a  coherent  combining  architecture 
can  be  used  for  DE  applications.  These  issues  include  the  consideration  of  sources  with  a  finite 
spectral  line  width  and  the  dephasing  effect  of  atmospheric  turbulence. 

Coherent  combining  has  proven  to  be  effective  in  situations  with  very  low-power  lasers  and 
weak  turbulence  [4] .  However,  these  conditions  are  not  applicable  for  DE  systems  which  require 
high-power  lasers  and  must  be  effective  in  conditions  of  moderate  and  strong  atmospheric 
turbulence.  High-power  fiber  lasers  have  significantly  broader  linewidths  due  to  stimulated 
Brillouin  scattering,  Doppler  shift,  self-phase  modulation,  and  Raman  broadening  in  the  gain 
medium.  Stimulated  Brillouin  scattering  (SBS)  is  the  main  contributor  to  the  large  linewidths  in 
high-power  fiber  lasers.  The  SBS  instability  results  in  a  threshold  power  level  for  fiber  and  slab 
lasers.  To  increase  the  threshold  power  level,  the  line  width  is  intentionally  broadened. 

The  full  width  half  max  (EWHM)  of  the  linewidth  for  high-power  fiber  lasers  (multi-kW)  is 
typically  Av  ~  lOGHz  and  the  coherence  time  is  tc  □  1/Al>  =  0.1  nsec  (for  all  single  peaked 
lineshapes)  [9].  Phase  locking  of  lasers  requires  measuring  the  output  phase  information  and 
applying  the  corrective  phase  to  the  individual  lasers.  To  be  effective,  the  time  scale  required  for 
this  process  should  be  less  than  the  laser  coherence  time.  The  coherence  time  is  the  characteristic 
time  over  which  the  phase  and  intensity  randomly  vary  and  is  due  solely  to  the  finite  spectral 
linewidth  in  the  gain  medium  and  the  statistical  nature  of  the  emission  from  the  atoms  or 
molecules  [8].  Stated  simply,  if  the  process  of  measuring  the  output  phase  and  applying  the 
corrective  phase  takes  longer  than  the  coherence  time,  then  the  output  phase  has  changed  before 
the  corresponding  corrective  phase  can  be  applied.  To  the  best  of  our  knowledge,  there  do  not 
currently  exist  instruments  that  can  operate  at  rates  comparable  to  the  EWHM  of  the  linewidth 
for  high-power  fiber  lasers.  However,  in  the  following  work  we  assume  that  there  exist 
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instruments  that  can  control  the  phase  over  time  scales  shorter  than  t^Ulj  ls.v  =  0.1  nsec . 
Furthermore,  we  assume  that  these  instruments  operate  with  zero  error. 

As  previously  stated,  the  power  spectral  linewidth  associated  with  the  gain  medium  results  in 
temporal  fluctuations  of  the  intensity  and  phase  at  the  output.  The  nature  of  the  temporal 
fluctuations  depends  on  the  shape  of  the  power  spectral  linewidth.  As  shown  in  Appendix  B,  the 
fluctuations  associated  with  a  Lorentzian  power  spectral  density  are  characteristically  different 
than  the  fluctuations  associated  with  a  Gaussian  power  spectral  density  due  to  the  extended 
wings  of  the  Lorentzian  spectrum.  A  Lorentzian  spectrum  can  result  from  collisional  broadening, 
while  a  Gaussian  spectrum  can  result  from  Doppler  broadening  or  local  random  atom 
environments  in  a  solid.  Regardless  of  the  lineshape,  the  standard  deviation  of  intensity  is 
approximately  equal  to  the  average  intensity.  In  addition,  the  output  phase  of  the  field  randomly 
fluctuates  over  2n  radians  on  time  scales  down  to  the  coherence  time.  Experimental 
observations  on  the  lineshapes  of  high-power  fiber  lasers  have  revealed  lineshapes  that  are 
approximately  Lorentzian. 

In  this  paper  we  discuss  the  physical  processes  associated  with  the  combining  and  propagation 
of  high -power  laser  beams.  We  compare  the  energy  delivered  to  a  target  for  the  case  of 
coherently  combined  and  incoherently  combined  laser  beams.  Through  simulation,  we 
demonstrate  that  in  vacuum  and  weak  turbulence,  the  effectiveness  of  coherent  combining  is 
limited  by  the  phase  fluctuations  occurring  due  to  the  broad  linewidth  of  high-power  lasers 
(Figures  2  and  3).  Furthermore,  the  advantages  of  coherent  combining  begin  to  diminish  in 
moderate  turbulence,  as  the  turbulence  induced  phase  distortions  become  dominant  (Figures  5 
and  6).  In  strong  turbulence  and  multi-km  propagation  ranges,  we  find  negligible  differences  in 
the  energy  on  the  target  between  coherently  and  incoherently  combined  laser  beams  (Figures  7 
and  8).  The  incoherently  combined  architecture  is  far  simpler  to  implement  and  both  beam 
combining  architectures  have  the  capability  of  employing  adaptive  optics  to  extend  the  range. 

We  begin  in  Section  II  by  describing  the  model  for  the  beam  director.  We  present  our 
propagation  simulation  results  in  Section  III  for  weak,  moderate,  and  strong  turbulence.  A 
discussion  of  our  findings  concludes  this  paper  in  Section  IV. 

II.  Beam  Director  Model 
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We  model  the  beam  direetor  as  a  set  of  nine  square  beams  which  we  refer  to  as  tiles.  The  side 
length  of  each  tile  is  a  and  the  distance  between  the  centers  of  adjacent  tiles  is  b  where  b>a.  For 
example,  a  beam  director  with  the  parameters  a=b  would  have  no  filling  factor  and  appear  as  a 
single  tile  with  side  length  3a  and  centered  at  (;t  =  0,  y  =  0) .  The  tile  arrangement  in  the  xy  plane  is 
shown  in  Fig.  la.  Tiles  are  indexed  by  the  indices  I  and  m  corresponding  to  the  x  and  y 
dimensions  respectively.  Indices  I  and  m  are  integers  defined  on  the  set  {-1,0,1} .  For  example,  the 
indices  (l  =  -\,m  =  -  1)  ,(/  =  0,  m  =  O) ,  and  (/  =  1,  m  =  l)  correspond  to  the  bottom  left,  center,  and  top 
right  tiles  respectively. 

Figure  lb  shows  a  side  view  of  the  beam  director,  i.e.  yz  plane.  Each  tile  is  positioned  with 
tip-tilt  correction  so  that  the  center  of  the  tile  propagates  to  the  point  [x  =  0,y  =  0,z  =  L)  in  vacuum 

where  L  is  the  distance  from  transmitter  to  the  target  and  L  □  b.  The  tilt  applied  to  a  tile  with 
index  m  is  defined  as  cpy  =  mbj L .  Similarly  the  tip  applied  to  a  tile  with  index  I  is  defined  as 

(p,=lblL. 


Figure  1:  Beam  director  geometry.  We  define  the  side  length  of  each  tile  as  a  and  the  distance 
between  the  centers  of  adjacent  tiles  as  b  where  b>a.  Each  tile  is  indexed  by  the  indices  I  and  m 
defined  on  the  set  {-1,0,1} . 

We  follow  the  method  presented  in  [10]  and  [11]  for  generation  of  a  laser  beam  with  arbitrary 
spectral  power  density.  The  method  entails  modeling  the  beam  as  radiation  from  a  large  number 
of  independent  radiators  that  radiate  at  discrete  frequencies.  We  consider  radiation  from  a  total  of 
N  frequencies  and  write  the  n*  frequency  as  +  nSco ,  where  o,,  is  the  center  frequency,  Sco 
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is  a  small  change  in  frequency,  and  n  is  an  integer  value  ranging  from  -Njl  to  Njl-l.  The  time 
dependence  of  the  electric  field  for  a  single  frequency  can  be  written  as 

=  +C-C-,  (1) 

where  and  are  random  variables.  As  a  consequence  of  the  central  limit  theorem,  and 
are  Gaussian  random  variables  with  zero  mean  and  variance  proportional  to  the  number  of 
radiators  at  frequency  co^ ,  which  is  determined  by  the  power  spectrum  of  radiation.  The 

implications  of  Gaussian  and  Lorentzian  power  spectra  are  discussed  in  Appendices  A  and  B. 

The  time  dependence  of  the  total  electric  field  from  all  frequency  components  is 

£(')  =  (')  =  iZ<“.  +»„)<■"'"■'  +C-C-  ■  (2) 

n  n 

Setting  =  OTq  +  nSco ,  the  total  electric  field  is  given  by 

E{t^  =  \  exp  ^  (a„  +  ib^ )  exp \-in5cot^  +  c.c. .  (3) 

n 

The  term  is  a  time  dependent  complex  value.  We  denote  this  value  as 

n 

w/l-l 

a{t^  +  ip{t^=  ^  +i/?„)exp[-OT<5'(a)t].  (4) 

n=-W/2 

Numerically,  t  is  a  sequence  of  discrete  values  such  that  t  —  mSt  where  m  is  an  integer  value 
ranging  from  0  to  A-1  and  St  is  to  be  determined. 

Njl-l 

a[mSt)  +  ip{^mSt)=  ^  (^a^+ip)e.xp[-inmScoSty  (5) 

n=~NI2 

If  we  choose  St  =  2nj {^NSco) ,  we  get  the  exact  definition  of  a  discrete  Fourier  transform  (DFT) 

^72-1 

a[mSt)  +  ip[mSt)=  ^  [a^+ip)exp[-i2n:nm/N^  (6) 

n=-W/2 

Thus  it  is  convenient  to  solve  for  (z(t)  and  with  the  use  of  a  DFT  on  the  sequence 

+  ib^ )  where  a  and  p  are  the  real  imaginary  part  of  the  resulting  DFT  respectively.  The  total 
electric  field  can  be  written  as  an  amplitude  and  phase  modulated  monochromatic  source 

£■(?)  =  i(a(t)  +  i/?(t))exp[-i0ot]  +  c.c.  .  (7) 
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Combining  the  time  dependence  described  by  Eq.  (7)  with  the  spatial  geometry  of  the  beam 
director,  we  write  the  electric  field  of  a  single  tile  as 

(x,  y,  z  =  0,  t)  =  i («,  „  (t)  +  (t)) rect  (x  - lb)xQci {y-mb) 

(8) 

where  rect  (4")  is  the  rectangular  function  defined  as 

0  if  1^1  >  a/2 

rect(^)  =  1/2  if  1^1  =  a/2  (9) 

1  if  1^1  <  a/2 

The  electric  field  of  the  beam  director  is  obtained  by  summing  over  all  tiles: 

/=1  m=l 

E(x,y,z  =  0,i)  =  ^  (0  +  *Am (0)rect(x-//7)rect(y-ml7) 

/=-lm=-l  .  (10) 

In  theory,  two  coherently  combined  tiles  would  have  a  phase  difference  of  zero  at  any 
instance  in  time.  However,  in  practice  this  is  generally  not  the  case  for  high-power  lasers. 

Instead,  coherent  combining  systems  result  in  a  reduction  of  the  root  mean  square  (RMS)  phase 
difference  between  tiles.  For  reference,  the  RMS  phase  difference  between  two  incoherently 
combined  tiles  is  2;r/^/i2  as  shown  in  Appendix  C.  To  simulate  a  situation  where  the  tiles  are 
coherently  combined  in  the  transmitter  plane,  we  generate  a  set  of  a,  „  (f)  and  „  (f)  in  which  the 
RMS  phase  difference  between  any  two  tiles  can  be  precisely  controlled.  This  procedure  for 
generating  a  set  of  correlated  a,  „  (t)  and  „  (t)  is  described  in  Appendix  C. 

III.  Propagation  in  Turbulent  Atmosphere 

We  simulate  the  propagation  of  the  tile  arrangement  described  in  Section  II  through  the 
atmosphere  for  coherently  and  incoherently  combined  tiles.  For  comparison,  we  include  the  case 
of  monochromatic,  phase  matched  tiles.  Following  our  previous  work  [12],  we  propagate  the 
beam  by  numerically  solving  the  paraxial  wave  equation.  Atmospheric  turbulence  is  modeled  as 
phase  screens  located  at  discrete  locations  along  the  z-axis.  For  details  on  the  split-step  phase 
screen  simulation  see  [12]. 
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We  use  parameters  that  are  typical  for  DE  applications.  Specifically,  the  center  wavelength  of 
radiation  in  vacuum  A,,  =  l|am ,  distance  to  the  target  L  =  5  km ,  tile  side  length  a  =  2  cm ,  distance 
between  tile  centers  b  =  2.2  cm ,  bandwidth  Av  =  A(y/2;r  =  37.5  GHz ,  and  RMS  phase  difference 
between  any  two  coherently  combined  tiles  =  2;r/6 .  We  consider  the  lineshape  to  be 
Lorentzian.  These  parameters  give  a  fractional  bandwidth  of  Aco/cog  □  1.3x10^  and  a  coherence 
time  of  □  8.5x10 sec .  We  consider  three  levels  of  atmospheric  turbulence,  =  0 , 

=10 and  =10 corresponding  to  vacuum,  moderate  turbulence  and  strong 
turbulence  respectively.  Using  the  Rytov  variance,  a\  =  10.5C^^^^^®L“''^ ,  as  a  measure  of  turbulence 
strength,  these  conditions  correspond  to  cr^  =  0 ,  g\=  6.3 ,  and  cr^  =  63 .  We  use  the  Fried 
parameter,  To  =  .184(c„^l/^^)  as  a  measure  of  the  transverse  coherence  length.  Strictly  speaking, 
Yg  corresponds  to  the  diameter  of  a  circular  area  in  which  the  RMS  phase  of  a  plane  wave  is  1 
radian, 
a)  Vacuum 

To  observe  effects  exclusively  associated  with  the  linewidth  of  the  laser  sources,  we  begin  by 
considering  propagation  in  vacuum  corresponding  to  =  0 ,  cr^  =  0 ,  and  ^^=00.  Figure  2 
displays  the  intensity  profiles  at  z  =  5  km  for  the  incoherently  combined  beam  (a),  coherently 
combined  beam  (b),  and  monochromatic,  phase  matched  beam  (c).  The  intensity  profiles  are 
averaged  over  a  time  scale  much  greater  than  the  coherence  time,  specifically  6  ns.  The 
incoherently  combined  tiles  diffract  independently  resulting  in  a  radially  symmetric, 
approximately  Gaussian  intensity  profile.  The  energy  from  the  incoherently  combined  beam  is 
spread  over  a  much  larger  area  than  for  the  coherently  combined  beams.  The  coherently 
combined  beam  produces  a  far  field  intensity  pattern  similar  to  the  monochromatic  beam. 
However,  there  are  considerable  differences.  The  coherently  combined  beam  has  diffracted  more 
than  the  monochromatic  beam  which  results  in  more  beam  spreading  and  a  lower  on  axis 
intensity. 
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Figure  2:  Intensity  profiles  for  propagation  through  vacuum  at  z  =  5  kmfor  the  incoherently 
combined  beam  (a),  coherently  combined  beam  (b),  and  phase  matched  monochromatic  beam 
(c).  Each  tile  of  the  incoherently  combined  beam  diffracts  independently  resulting  in  significant 
spreading.  The  coherently  combined  beam  has  an  initial  RMS  phase  difference  =  2;t/6 

between  any  two  tiles.  The  far- field  intensity  profile  of  the  coherently  combined  beam  resembles 
that  of  the  monochromatic  beam. 


Figure  3  displays  the  intensity  profile  along  the  transverse  x-axis  for  the  incoherent  beam 
(blue),  coherent  beam  (green),  and  monochromatic,  phase  matched  beam  (red).  Intensity  is 
normalized  by  the  on-axis  intensity  of  the  monochromatic,  phase  matched  beam.  The  on  axis 
intensity  of  the  coherently  combined  beam  is  approximately  6.4  times  greater  than  the  on-axis 
intensity  of  the  incoherently  combined  beam.  However,  the  decrease  in  on-axis  intensity  of  the 
coherently  combined  beam  in  comparison  to  the  monochromatic  beam  is  a  fundamental 
limitation  caused  by  the  random  temporal  characteristics  of  the  phase  associated  with  each 
individual  tile. 
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Figure  3:  Vacuum  intensity  profiles  along  x-axis  at  z  =  5  km  for  the  incoherently  combined  beam 
(blue),  coherently  combined  beam  (green),  and  phase  matched  monochromatic  beam  (red). 
Intensity  is  normalized  by  the  on-axis  intensity  of  the  monochromatic  beam.  The  on  axis 
intensity  of  the  incoherently  and  coherently  combined  beams  are  approximately  11%  and  71%  of 
the  on-axis  intensity  of  the  monochromatic  beam  respectively. 

A  common  metric  used  in  DE  applications  is  power  in  the  bucket  (PIB),  which  describes  the 
amount  of  power  contained  in  a  specific  area.  We  calculate  the  PIB  centered  at  the  origin  as  a 
function  of  bucket  radius  for  the  time  averaged  intensity  profiles.  The  results  are  displayed  in 
Fig.  4  where  the  PIB  is  normalized  by  the  total  power.  Again,  incoherent,  coherent,  and 
monochromatic  results  are  denoted  by  blue,  green,  and  red  respectively.  The  coherently 
combined  beam  and  monochromatic  beam  each  have  a  distinct  central  lobe  with  a  radius  of 
approximately  7.5  cm.  At  a  radius  of  7.5  cm,  the  incoherent  and  coherently  combined  beams 
contain  37%  and  78%,  respectively,  of  the  power  contained  by  the  monochromatic  beam.  Unlike 
the  coherently  combined  and  monochromatic  beam,  the  incoherently  combined  beam  remains 
radially  symmetric  at  radii  greater  than  the  central  lobe  radius.  Consequently,  the  PIB  of  the 
incoherently  combined  beam  increases  more  rapidly  than  coherently  combined  beam  after  a 
radius  of  7.5  cm.  At  a  radius  of  20  cm,  the  contained  power  is  approximately  equal  for  all  beams. 
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Figure  4:  :  Power  in  the  bucket  (PIB)  at  z  =  5  km  through  vacuum  as  a  function  of  bucket  radius 
for  the  incoherently  combined  beam  (blue),  coherently  combined  beam  (green),  and  phase 
matched  monochromatic  beam  (red).  The  coherently  combined  beam  and  monochromatic  beam 
contain  significantly  more  power  near  the  axis  due  to  the  distinct  central  lobe.  The  contained 
power  within  a  radius  of  20  cm  is  approximately  equal  for  all  beams. 


b)  Moderate  turbulence 

Now  we  consider  propagation  through  moderate  turbulence  corresponding  to 
=  10 ,  al  =  6.3 ,  and  =1.76  cm .  All  Other  parameters  remain  the  same.  An  ensemble 
average  is  performed  over  the  intensity  profiles  from  100  independent  instances  of  turbulence. 
The  intensity  profile  for  each  instance  of  turbulence  is  averaged  over  many  coherence  time 
intervals,  6  ns,  to  account  for  the  fluctuating  intensity  and  phase.  The  average  intensity  profiles  at 
z  =  5  km  are  displayed  in  Fig.  5  for  the  incoherently  combined  beam  (a),  coherently  combined 
beam  (b),  and  phase  matched,  monochromatic  beam  (c).  Through  comparison  of  the  average 
intensity  profiles  we  can  observe  a  slight  advantage  of  coherent  beam  combining  over  incoherent 
beam  combining.  These  advantages  include  a  smaller  spot  size  and  higher  maximum  intensity. 
However,  the  effects  of  coherent  beam  combining  are  nearly  insignificant  when  compared  to  the 
impact  of  moderate  turbulence  as  shown  in  Fig.  6. 
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Figure  5:  Intensity  profiles  for  propagation  through  moderate  turbulence,  =  10  m  ,  at 
z  =  5  kmfor  the  incoherently  combined  beam  (a),  a  coherently  combined  beam  (b),  and  phase 
matched  monochromatic  beam  (c).  The  coherently  combined  beam  and  monochromatic  beam 
have  smaller  spot  size  and  higher  on-axis  intensity  than  the  incoherently  combined  beam. 

Figure  6  displays  the  intensity  profile  along  the  x-axis  for  the  incoherent  (blue),  coherent 
(green),  and  monochromatic  (red)  beam  profiles.  Again,  the  intensity  is  normalized  by  the  on- 
axis  intensity  of  the  monochromatic  beam  in  vacuum  (black).  The  on-axis  intensity  of  the 
coherently  combined  beam  in  moderate  turbulence  is  approximately  25%  of  the  on-axis  intensity 
of  the  coherently  combined  beam  in  vacuum.  Turbulence  has  a  profound  impact  on  the 
coherently  combined  beam  in  comparison  to  the  incoherently  combined  beam.  The  on-axis 
intensity  of  the  incoherently  combined  beam  in  moderate  turbulence  is  approximately  78%  of  the 
on-axis  intensity  of  the  incoherently  combined  beam  in  vacuum. 
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Figure  6:  Average  intensity  profiles  in  moderate  turbulence,  =  10  m  ,  along  x-axis  at 
z  =  5  kmfor  the  incoherently  combined  beam  (blue),  coherently  combined  beam  (green),  and 
phase  matched  monochromatic  beam  (red).  Intensity  is  normalized  by  the  on-axis  intensity  of  a 
monochromatic  beam  propagated  through  vacuum.  Propagation  through  vacuum  is  denoted  by 
the  black  curve.  The  on-axis  intensity  of  the  coherently  combined  is  approximately  twice  the  on- 
axis  intensity  of  the  incoherently  combined  beam. 

We  also  find  that  moderate  turbulence  significantly  impacts  the  advantage  of  PIB  for 
coherently  combined  beam.  Figure  7  displays  the  PIB  as  a  function  of  radius  at  z  =  5  kmfor 
beams  propagating  through  moderate  turbulence.  The  coherently  combined  beam  and 
monochromatic  beam  no  longer  have  a  distinct  central  lobe  as  in  the  case  of  propagation  through 
vacuum.  As  a  result,  the  PIB  of  each  beam  has  a  similar  trend.  The  coherently  combined  beam 
and  monochromatic  beam  still  contain  more  power  near  the  axis  than  the  incoherently  combined 
beam,  but  the  difference  is  not  as  significant.  The  advantage  of  coherent  combining  appears 
minor  when  compared  to  the  PIB  to  a  phase  matched,  monochromatic  beam  in  vacuum.  Even  in 
moderate  turbulence,  it  is  clear  that  the  phase  distortions  caused  by  atmospheric  turbulence  are 
dominant  over  the  effects  of  coherent  beam  combining. 
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Figure  7:  Power  in  the  bucket  (PIB)  at  z  =  5  km  through  moderate  turbulence,  Cl  =10  m  ,  as  a 
function  of  bucket  radius  for  the  incoherently  combined  beam  (blue),  coherently  combined  beam 
(green),  and  phase  matched  monochromatic  beam  (red).  For  reference,  the  PIB  of  the 
monochromatic,  phase  matched  beam  propagated  through  vacuum  is  denoted  by  the  dashed 
black  curve.  The  coherently  combined  beam  and  monochromatic  beam  contain  more  power  near 
the  axis  than  the  incoherently  combined  beam,  however  the  difference  is  not  as  significant  as  in 
vacuum. 


c)  Strong  turbulence 

Finally,  we  consider  propagation  through  strong  turbulence  corresponding  to  Cl  =  10 , 
al=63,  and  =  .44  cm .  All  Other  parameters  remain  the  same.  Again,  an  ensemble  average  is 
performed  over  the  intensity  profiles  associated  with  100  independent  instances  of  turbulence 
and  the  intensity  profile  for  each  instance  of  turbulence  is  averaged  over  6  ns.  The  average 
intensity  profiles  at  z  =  5  km  are  displayed  in  Fig.  8  for  the  incoherently  combined  beam  (a), 
coherently  combined  beam  (b),  and  phase  matched,  monochromatic  beam  (c).  All  three  intensity 
profiles  are  nearly  indistinguishable  from  each  other.  The  intensity  profiles  of  the  coherently 
combined  beam  and  monochromatic  beam  have  no  resemblance  to  the  vacuum  far- field  pattern 
displayed  in  Fig.  2b,c. 
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Figure  8:  Intensity  profiles  for  propagation  through  strong  turbulence,  C]  =  10  m  at  z  =  5  km 
for  the  incoherently  combined  beam  (a),  coherently  combined  beam  (b),  and  phase  matched 
monochromatic  beam  (c).  The  intensity  profiles  are  nearly  indistinguishable. 


Figure  9  displays  the  intensity  profile  along  the  x-axis  for  the  incoherent  beam  (blue), 
coherent  beam  (green),  and  monochromatic  beam  (red).  Again,  the  intensity  is  normalized  by  the 
on-axis  intensity  of  the  monochromatic  beam  in  vacuum.  Atmospheric  turbulence  is  dominant 
over  the  method  of  beam  combining;  consequently,  the  shape  of  each  intensity  profile  is  nearly 
identical.  The  main  noticeable  difference  is  that  the  coherently  combined  beam  and 
monochromatic  beam  have  slightly  higher  spatial  fluctuations  in  the  time  averaged  intensity. 
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Figure  9:  Intensity  profiles  in  strong  turbulence,  =  10  m  ,  along  x-axis  at  z  =  5  kmfor  the 
incoherently  combined  beam  (blue),  coherently  combined  beam  (green),  monochromatic  beam 
(red).  Intensity  is  normalized  by  the  on-axis  intensity  of  a  monochromatic  beam  in  vacuum.  The 
intensity  profiles  are  nearly  indistinguishable. 


Figure  10  displays  the  PIB  as  a  function  of  radius  for  the  incoherent  beam  (blue),  coherent 
beam  (green),  and  monochromatic  beam  (red).  The  initial  condition  of  the  phase  at  z  =  0  has  very 
little  discernible  impact  on  the  PIB  at  z  =  5  km  in  strong  turbulence.  The  advantages  of  coherent 
beam  combining  that  were  apparent  in  vacuum  have  completely  diminished  in  strong  turbulence. 
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Figure  10:  Power  in  the  bucket  (PIB)  at  z  =  5  km  through  strong  turbulence,  Cl  =  10  m  ,  as  a 
function  of  bucket  radius  for  the  incoherently  combined  beam  (blue),  coherently  combined  beam 
(green),  and  phase  matched  monochromatic  beam  (red).  All  beams  deliver  approximately  the 
same  power. 


IV.  Discussion 

We  have  demonstrated  that  the  linewidth  of  the  source  places  a  fundamental  limit  on  the 
ability  to  coherently  combine  beams.  The  large  linewidth  associated  with  high-power  lasers  must 
be  carefully  considered  if  they  are  to  be  coherently  combined.  However,  coherently  combining 
beams  at  the  transmitter  plane  has  no  benefit  for  DE  systems  in  conditions  of  strong  turbulence. 
We  have  shown  that  the  initial  RMS  phase  difference  between  tiles  has  negligible  impact  on  the 
quality  of  the  beam  after  propagation  through  multiple  kilometers  of  moderate  or  strong 
atmospheric  turbulence.  In  strong  turbulence,  coherent  beam  combining  at  the  transmitter  plane 
merely  adds  to  the  complexity  of  the  transmitter.  Situations  of  strong  turbulence  and  multi -km 
distances  require  adaptive  optics  to  compensate  for  the  phase  distortions  caused  by  atmospheric 
turbulence.  Depending  on  the  implementation  of  the  adaptive  optics,  the  linewidth  of  high-power 
lasers  may  need  to  be  taken  into  consideration.  For  instance,  coherent  combining  at  the  target  can 
be  considered  an  adaptive  optics  method  because  it  partially  compensates  for  distortions  caused 
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by  atmospheric  turbulence.  In  this  case,  the  line  width  of  high-power  lasers  places  even  greater 
limitations  on  coherent  combining  due  to  the  transit  time  of  light  to  the  target  and  back  to  the 
receiver.  If  the  transit  time  is  longer  than  the  coherence  time  of  the  lasers,  then  it  is  inconceivable 
to  coherently  combine  the  lasers  at  the  target  without  first  phase  matching  the  beams  at  the 
transmitter. 

There  are  multiple  limitations  that  the  linewidth  of  a  laser  places  on  the  ability  to  coherently 
combine  beams.  Due  to  the  dominant  effect  of  atmospheric  turbulence,  it  is  not  effective  to 
coherently  combine  lasers  at  the  transmitter  plane.  Incoherent  combining  of  lasers  is  a  much 
simpler  approach  with  comparable  results  in  moderate  to  strong  turbulence.  A  more  effective 
approach  for  delivering  energy  to  a  target  in  moderate  to  strong  turbulence  may  be  in  using 
adaptive  optics  solutions  to  compensate  for  the  turbulent  distortions. 
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Appendix  A 

For  a  continuous  power  spectrum  of  radiation,  {co) ,  the  variance  of  and  in  Eq.  (1)  are 
Var(a„ )  =  Var(Z7„ )  =  Scoa^  (o„ ) .  Here  we  consider  a  Gaussian  power  spectrum,  ,  and  a  Lorentzian 
power  spectrum,  al ,  defined  as 


{(o)  = 


exp 


2r? 


2  \ 


(Al) 


al  (co)  =  - 
n 


(A2) 


respectively,  where  and  are  parameters  that  determine  the  width  of  each  power  spectrum. 
The  power  spectra  are  normalized  such  that 

+O0  +O0 

I  cTq  {co)dco  =  I  cTl  {G))dco  =  1 .  (A3) 


The  line  width,  Azy ,  is  defined  as  the  full  width  at  half  maximum  (FWHM)  of  the  power 
spectrum.  The  FWHM  of  a  Gaussian  function  is  defined  as  =  2(21n2)’^^  and  the  FWHM 
of  a  Forentzian  function  is  defined  as  Azy^  =  2r^ .  For  comparison  of  the  two  spectra,  it  is 
convenient  to  pick  r^,  and  such  that  Aco^  =  Aco^  =  Aco .  For  example,  =  lOOGHz  and 
=  (21n 2)''^^  Fq  =  1 17.7GHz  which  gives  the  relation  Azy^  =  Azy^  =Aco  =  235.5GHz .  The  coherence 
time  is  defined  as  =  2/ Aco .  The  linewidth  ratio  for  this  example  is  Aco/ =  1.24x  10^  .  Figure 
Al  displays  a  comparison  of  two  power  spectra  for  these  parameters. 
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Figure  Al:  Comparison  of  a  Gaussian  and  Lorentzian  power  spectrum  when  the  line  widths  are 
matched,  i.e.  =  A® .  The  Gaussian  power  spectrum  contains  98%  of  the  total  power 

within  Afij  of  .  The  Lorentzian  power  spectrum  contains  70%  of  the  total  power  within  A«j  of 

%■ 

The  Lorentzian  power  spectrum  decays  more  slowly  than  the  Gaussian  power  spectrum  and 
thus  a  greater  fraction  of  the  total  power  is  contained  at  frequencies  far  from  co^ .  For  example, 
the  Gaussian  power  spectrum  contains  98%  of  the  total  power  within  hco  of  ,  while  the 
Lorentzian  power  spectrum  contains  70%  of  the  total  power  within  Ac;  of  .  Appendix  B 
illustrates  this  importance  by  comparing  the  random  fluctuations  in  the  intensity  and  phase  for 
Gaussian  and  Lorentzian  power  spectral  densities. 

Appendix  B 

Here  we  illustrate  the  random  fluctuations  in  the  intensity  and  phase.  We  examine  the 
behavior  of  the  intensity  and  phase  for  a  single  tile  at  z  =  0  and  on  time  scales  comparable  to  the 
coherence  time,  C  •  The  intensity  of  a  tile  is  defined  as 

I(l)  =  ^{E{tf)  =^{a‘{t)+fi^t)),  (Bl) 

4;r  '  ! t=ini(o^  ' 
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where  the  angled  brackets  indicate  a  temporal  average  performed  over  a  time  scale  of  2;r/®o  .  It 
is  convenient  to  define  a  normalized  intensity  as  /„  (?)  =  (?)  +  (?) .  Due  to  the  statistical 

properties  of  a  and  p  we  have  the  relation  ^  ^  phase  of  a  tile  is  defined  as 

(^(?)  =  arctan[/?(t)/(z(t)] .  (B2) 

We  note  that  strictly  speaking,  the  arctan  function  has  a  range  of  (-;t/2,;t/2)  .  Numerically,  we 
utilize  the  atan2  function  so  that  the  phase  of  a  tile  is  defined  on  the  interval  (-;r,;r] . 

We  begin  by  considering  a  tile  with  a  Gaussian  power  spectrum  defined  by  Eq.  (Al).  Figure 
B 1  shows  the  intensity  fluctuations  at  z  =  0  for  a  single  tile  with  a  Gaussian  power  spectrum.  The 
intensity  fluctuates  randomly,  but  still  appears  as  a  smooth  function  in  time.  Major  fluctuations 
occur  over  intervals  greater  than  the  coherence  time  as  displayed  in  figure  Fig.  Bl.  Calculating 
the  long  term  average  intensity  for  this  specific  instance  of  a  and  p  gives  the  result 


Figure  Bl:  Intensity  fluctuations  for  a  single  tile  with  Gaussian  power  spectrum.  The  intensity 
has  random  fluctuations  when  observed  on  time  scales  comparable  to  the  coherence  time.  The 
standard  deviation  of  intensity  is  equal  to  the  average  intensity. 

Figure  B2  shows  the  phase  fluctuations  at  z  =  0  for  a  single  tile  with  a  Gaussian  power  spectrum. 
Similar  to  the  intensity,  the  phase  fluctuates  randomly  when  observed  on  a  time  scale 
comparable  to  the  coherence  time. 


21 


3 


1 


-31-  •  ,1  _  'si,  ,  1 _ 

0  10  20  30  40  50  60  70  80  90  100 

t  / [arb.  units] 

Figure  B2:  Phase  fluctuations  for  a  single  tile  with  Gaussian  power  speetrum.  The  phase  fluctuates 
randomly  when  observed  on  time  seales  eomparable  to  the  eoherenee  time. 

The  intensity  and  phase  as  a  function  of  time  is  characteristically  different  when  the  tile  has  a 
Lorentzian  power  speetrum  deseribed  by  Eq.  (A2).  Figure  B3  shows  the  intensity  fluctuations  at 
z  =  0  for  a  single  tile  with  a  Lorentzian  power  spectrum.  Similar  to  the  ease  of  the  Gaussian 
power  spectrum,  large  seale  fluetuations  oeeur  over  times  of  a  few  eoherenee  time  intervals. 
However  the  intensity  fluctuations  do  not  appear  as  a  smooth  funetion.  Regardless,  performing  a 
long  time  average  gives  the  result  □  -98 ,  which  is  close  to  the  expected  value  of  1. 


t  / [arb.  units] 


Figure  B3:  Intensity  fluctuations  for  a  single  tile  with  a  Lorentzian  power  speetrum.  The  standard 
deviation  of  intensity  is  equal  to  the  average  intensity. 
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Figure  B4  shows  the  phase  fluctuations  at  z  =  0  for  a  single  tile  with  a  Lorentzian  power 
spectrum.  The  phase  displays  the  same  random  fluctuations  over  time  scales  comparable  to  the 
coherence  time. 
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Figure  B4:  Phase  fluctuations  for  a  single  tile  with  a  Lorentzian  power  spectrum. 

Appendix  C 

Ideally,  the  phase  of  two  coherently  combined  beams  would  fluctuate  synchronously  and  have 
a  phase  difference  of  zero  at  any  instance  in  time.  However,  in  practice  this  is  generally  not  the 
case,  especially  for  high-power  lasers.  Instead,  coherent  combining  systems  result  in  a  reduction 
of  the  root  mean  square  (RMS)  phase  difference  between  beams.  To  simulate  a  situation  where 
the  tiles  are  coherently  combined  in  the  transmitter  plane,  we  develop  a  method  to  generate  a  set 
of  a,  „  (f)  and  (?)  in  which  the  RMS  phase  difference  between  any  two  tiles  can  be  controlled. 
As  an  illustration,  we  consider  2  tiles  whose  amplitude  and  phase  are  modulated  by  a,  (?) ,  /?,  (?) , 

(0 ,  and ydj  (?) .  These  variables  are  generated  through  the  Fourier  transform  of  the  Gaussian 
random  variables  where  we  have  added  the  superscript  to  distinguish  between 

tiles.  If  we  assume  that  perfect  coherent  combining  is  implemented  for  all  frequency  components 
within  a  certain  frequency  bandwidth,  |dy„  -  (Wo|  <  Q ,  where  Q  is  the  half  width  of  the  coherent 
combining  instrument  bandwidth,  then  a®  =  af  ’  and  Z?®  =  ’  for  \n\  <  Q/do .  In  other  words, 

random  variables  associated  with  frequency  components  contained  within  the  coherent 
combining  bandwidth  are  shared  between  the  tiles.  Conversely,  random  variables  associated  with 
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frequency  components  outside  of  the  coherent  combining  bandwidth  are  uncorrelated  between 
the  tiles.  The  coherent  combining  bandwidth  is  determined  by  the  time  scale  over  which  the 
coherent  combining  instruments  operate.  We  assume  that  the  instruments  have  zero  error  and 
thus  can  precisely  phase  match  all  frequency  components  within  the  instrument  bandwidth. 

We  define  ^,,^2  the  phase  of  tile  1,  phase  of  tile  2,  and  the  phase  difference 

between  the  two  tiles  respectively.  In  terms  of  (f) ,  (f) ,  (?) ,  and  (?)  we  have  the  relations 


(t)  =  arctan[/?j(t)/aj(t)]  ,  (Cl) 

^2(t)  =  arctan[/?2(t)/a2(t)]  ’ 

=  ■  (C3) 

The  RMS  phase  difference  between  the  two  tiles  is  defined  as: 

( 1 

AA«=  T7Ea«>’(')  ■  (C4) 

t  y 


We  can  analytically  calculate  for  the  cases  of  incoherently  combined  tiles  and  perfectly 
phase  matched  tiles.  If  the  two  tiles  are  incoherently  combined,  then  A^is  a  random  variable 
with  a  uniform  distribution  on  the  interval  [0,;r]  and  .  From  the  properties  of  a 

uniform  distribution  we  know  the  expected  value,  e[a^]  =  njl ,  and  variance,  Var[A^]  =  jn. 
Using  the  definition  of  variance,  Var[A']  =  E[A'^]-E[A']Y  we  calculate  E  /3and 

.  On  the  other  hand,  if  the  two  tiles  are  perfectly  phase  matched  then  is  uniquely 
zero  and  clearly  =  0  .  These  simple  situations  of  incoherent  combining  and  perfect  phase 
matching  correspond  to  Q  =  0  and  Q  =  co  respectively.  Coherently  combining  monochromatic 
sources  is  an  example  of  perfect  phase  matching.  We  calculate  intermediary  values  A^^^  (Q)  for 

both  Gaussian  and  Lorentzian  power  spectra  displayed  in  Fig.  Cl. 

For  a  Gaussian  power  spectrum,  the  RMS  phase  difference  between  the  two  tiles  is  effectively 
zero  when  frequency  components  within  Q  =  2Att>  of  are  matched.  However,  the  result  differs 
greatly  when  the  two  tiles  have  a  Lorentzian  power  spectrum.  Since  a  higher  percentage  of 
power  is  contained  at  frequencies  far  from  ,  A^^„^  (q)  decays  much  slower.  From  Fig.  Cl,  it  is 

clear  that  it  becomes  very  difficult  to  achieve  a  small  RMS  phase  difference  between  two  tiles 
with  a  Lorentzian  power  spectrum. 
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Figure  Cl:  The  RMS  phase  difference  between  two  tiles  with  a  Gaussian  power  spectrum  (blue) 
and  Lorentzian  power  spectrum  (red)  as  a  function  of  instrument  bandwidth. 

Recalling  that  Q  is  the  half  width  of  the  coherent  combining  angular  frequency  bandwidth, 
we  can  use  Fig.  Cl  to  calculate  the  rate  at  which  the  instruments  must  operate  to  achieve  any 
arbitrary  RMS  phase.  For  example,  to  achieve  coherent  beam  combining  with  an  RMS  phase 
difference  of  2;r/6  would  require  instruments  that  operate  at  a  rate  of  2(.48A(w)/2;t  »  36  GHz  and 
2Affl/2;?r  =  75  GHz  for  a  Gaussian  and  Lorentzian  linewidth  respectively.  In  the  following 
simulations,  we  consider  the  propagation  of  coherently  combined  beams  with  a  Lorentzian 
power  spectrum. 
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